"0","```r
# fit the models
scpd1.geom <- fitsad(scpdSM[1,scpdSM[1,]>0], 'geom')
scpd1.ls <- fitsad(scpdSM[1,scpdSM[1,]>0], 'ls')
scpd1.lnorm <- fitsad(scpdSM[1,scpdSM[1,]>0], 'lnorm')
# compare the fit of each model
AICtab(scpd1.geom, scpd1.ls, scpd1.lnorm, base=TRUE)
```"
"1","            AIC  dAIC df
scpd1.ls    68.6  0.0 1 
scpd1.lnorm 74.3  5.7 2 
scpd1.geom  79.4 10.8 1 
"
"0","```r
# transform to appropriate class for plotting
scpd1.geom.oc <- octavpred(scpd1.geom)
scpd1.ls.oc <- octavpred(scpd1.ls)
scpd1.lnorm.oc <- octavpred(scpd1.lnorm)
# plot the three models against one another and the data
plot(scpd1.oc)
lines(scpd1.geom.oc, col=\"red\")
```"
"0","```r
lines(scpd1.ls.oc, col=\"blue\")
lines(scpd1.lnorm.oc, col=\"green\")
```"
"0","```r
legend(\"topright\", c(\"geometric\", \"Logseries\", \"Lognormal\"), lty=1, col=c(\"red\",\"blue\", \"green\"))
```"
